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ABSTRACT 

A recently developed technique is presented for thermographic detection of flaws in composite materials by performing 
temperature measurements with fiber optic Bragg gratings, individual optical fibers with multiple Bragg gratings 
employed as surface temperature sensors were bonded to the surfaces of composites with subsurface defects. The 
investigated structures included a 10-ply composite specimen with subsurface delaminations of various sizes and depths. 
Both during and following the application of a thermal heat flux to the surface, the individual Bragg grating sensors 
measured the temporal and spatial temperature variations. The data obtained from grating sensors were analyzed with 
thermal modeling techniques of conventional thermography to reveal particular characteristics of the interested areas. 
Results were compared with the calculations using numerical simulation techniques. Methods and limitations for 
performing in-situ structural health monitoring are discussed. 
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1. INTRODUCTION 


Fiber optic sensors have been extensively studied and proposed for temperature and strain sensing in heath monitoring 
systems of aerospace structures and materials as well as many other applications 1 " 3 Compared to other sensors, Fiber 
optic sensors have the advantages of being lightweight and flexible, and requiring simpler wiring especially for 
distributed sensing. For an extensive heath monitoring system distributed sensors are an important requirement. Most 
efforts have focused on using distributed fiber optic strain sensing systems, with limited consideration given to 
temperature sensing systems. For fiber optic temperature sensing, the conventional techniques using Raman scattering 
can detect the average temperature over a long distance within a single fiber. Flowever the spatial resolutions for 
temperature readings-averaging over a length of fiber-are typically about one meter. Therefore they are not suitable for 
measuring the thermal responses at specific locations of the investigated materials. 

Another more attractive alternative is using fiber Bragg gratings (FBGs) for distributed temperature sensing (DTS). 
FBGs written in an optical fiber have a physical length of only a few millimeters. This size is nearly a point sensor. 
Flowever, if multiple FBGs are written at different (multiplexing) wavelengths and read with a time domain 
demodulation system, the total number of these FBGs in a singe fiber is limited due to the finite bandwidth of the laser in 
the detector system. To solve the above problems of spatial resolution and limited sensors, recently a new technique for 
fiber optic DTS has been developed using low reflectivity FBGs. Under this technique low reflectivity FBGs were 
written at the same nominal wavelength and read with a frequency domain demodulation system. 4 ' 5 This technique 
allows a single fiber to contain hundreds of FBGs, employed as temperature and strain sensors. 

An effective fiber optic DTS system has some advantages compared to a distributed strain sensing system for structure 
health monitoring. For a strain sensing system, in general, it is required to mechanically excite (or stress) the investigated 
structure or material. 6 For some cases this is impractical, if not impossible. Obviously, using a DTS system requires the 
investigated material to be thermally excited. Flowever, the temperature change of the material from the required 
thermal excitation can be only a few degrees or no more than tens of degrees. Such temperature variation is comparable 
with the temperature changes excited for conventional thermography techniques. 

Conventional thermographic techniques, in general, utilize a flash or quartz lamp as a heating source and an infrared 
imager to detect the thermal response of the investigated material. The heating pulse duration typically ranges from a 
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fraction of a second to a few seconds. The IR imager contains an array of several hundred by hundred detectors. These 
thermographic techniques are capable of large area inspection of aerospace structures and materials for their reliability 
and safety. One of the particularly important areas is the inspection of graphite fiber reinforced composite materials 
because they are being increasingly used as primary structures due to their high stiffness and strength to weight ratio. Of 
particular interest is the detection of delaminations that can appreciably reduce the compressive strength of a composite. 
Recently Winfree et al. successfully developed a physical model of two-layered systems for accurate reduction of the 
temporal thermal response of a composite with fabricated delaminations to the depths of the delaminations. 7 

In this paper we propose a new technique using a single optical fiber with multiple FBGs for the thermographic detection 
of flaws in materials and structures. The investigated structures included a 10-ply composite specimen with fabricated 
delaminations of various sizes and depths, similar to the one mentioned above. The optical fiber was bonded to the 
surface of the investigated composite. Both during and following the application of a thermal heat flux to the surface, 
the individual Bragg grating sensors measured the temporal and spatial temperature variations. The data obtained from 
individual FBGs were analyzed with thermal modeling to reveal particular characteristics within the area of interest. 


2. THEORY 


2.1 FBG as a temperature sensor 

In general, a fiber Bragg grating can be characterized by its Bragg wavelength, which is the center wavelength of the 
light reflected from the grating. The Bragg wavelength is given as 


~ 2n e ffA, ( 1 ) 

where n ett is the effective refractive index of the fiber core and A the grating period. For a fiber Bragg grating bonded 
onto or embedded in a polymeric substrate, a change in the temperature causes a change in the grating period due to not 
only the thermal expansion of the fiber but also the strain induced by thermal expansion of the substrate. In addition, the 
refractive index of the fiber core changes because of the thermo-optic effect. Combining all the above effects, the shift 
in the Bragg wavelength due to a temperature change, AT, is given as 

5?w3 / A. b = (1 — p e )S/ // + 5n c ff/n eff , (2) 

where p e is the photoelastic constant of the optical fiber and Sill is the thermally induced strain of the fiber. In general 
the photoelastic constant, the fractional index change, and the thermal expansion coefficients of the fiber and the 
substrate are temperature dependent and nonlinear, especially at low temperatures. However, for a finite temperature 
change, especially for the temperatures around and above the room temperature, Eq. (2) can be rewritten as a linear form 

as, 8 


5a b /Ab = [(1 -Pe)a p +q]AT, (3) 

where ^ = (l/n eff )(<5n eff 1ST), is the thermo-optic coefficient of the fiber and a p is the thermal expansion coefficient of 
the substrate. Eq. (3) takes into account the conditions when the thermal expansion coefficient and the physical 
dimension of the substrate are much greater than those of the optical fiber. This simple equation allows the FBG to 
perform as a temperature sensor of the substrate. 

2.2 Thermal response of one and two-layered systems 

It is possible to derive a one-dimension series solution for a single layer with periodic heating using the method of 
images, if both sides of the layer are thermally insulated. The series solution for the surface where the periodic heating is 
applied is 
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where the periodic heat flux is given by F cosfetf), K, k, and L are the thermal conductivity, diffusivity and thickness of 
the layer. There are two extremes of interest. First, for the case of a very thick layer or L» (k/2 /(o) la , the thermal 
response lags behind the periodic flux by 7t/4. Secondly, for a very thin layer or L« (k/2/ oj) 1 ' 1 , the thermal response lags 
behind excitation by n/2. All other thicknesses have phase differences that fall between these two limits. 

Real delaminations in composites are air gaps between two layers in the composite. It is well known that for long times, 
the response of the temperature over the delamination is not characteristic of a single layer response. The typical 
explanation for this effect is that lateral heat flow around the delamination diffuses the heat to the backside of the 
delamination. However if the thickness of the air gap is such that the heat flow across the gap is much larger than the 
heat flow around the delamination, the gap’s thermal resistance dominates the time response of the delamination. Note 
that if the thermal resistance of the delamination is sufficiently small, the time response at the front surface may not be 
significantly different than a single layer response. 

A simple analytic solution does not exist for the one-dimensional heat flow in a multilayered material. A solution does 
however exist in Laplace space for two layers of thickness It and l 2 coupled by an intermediate contact resistance (R). 
Since the configuration of interest is a composite with a delamination, the first and second layers are assumed to have the 
same thermal conductivity (K) and diffusivity (k). For the surface with the incident heating, the Laplace transform of the 
temperature impulse response is given by 7 
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where q is ^p/K and p is the coordinate in Laplace space. 

The flux input from the periodically cycled heat lamps is assumed to be 

/(0 = ./o(l-cos(®0), (6) 

where ft is the energy per area from the heating source and t 0 is the heating duration. While it is possible to include this 
into Eq. (5), the numerical inverse Laplace transform of the combined equation is found to be unstable for periodic heat 
fluxes at times greater than a 3 or 4 periods. Therefore to determine the thermal response for the periodic heat source, the 
inverse Laplace transform Eq. (5) is found numerically and this response is convolved with Eq. (6), 

T f (t) = f T ' (t - r)(l - cos(a>T))dr , (7) 
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where T‘(t ) is the inverse of Eq. (4). 


3. EXPERIMENT 

The low (smaller than a few tenths of one percent) reflectivity FBGs used in this research were written in situ, into the 
optical fiber drawn using the NASA Langley optical fiber draw tower. They were written with a pulsed KrF-excimer 
laser of 248 nm and a Talbot interferometer arrangement mounted on the tower. The interferometer consisted of a phase 
mask functioning as a beam splitter and a pair of mirrors used to recombine the split beams to form an interference 
pattern. An aperture was placed in the laser beam path to control the grating length and spatial profile. 9 FBG lengths for 
the present study were nominally 5 mm. The grating pitch written into the fiber could be adjusted by changing the 
relative angle of the two mirrors. These single-mode fibers were drawn from commercially available germanium-doped 
preforms of high numerical apertures. The drawn fibers with FBGs were coated with polyimide to thicknesses ranging 
from 11 to 16 micrometers, then ink-marked to show grating locations. 



These low reflectivity FBGs were interrogated using a frequency domain demodulation system shown in Fig. 1. In this 
system, the fiber coupler Ci and a pair of Faraday rotation mirrors (FRMs) form an in-fiber interferometer with an 
optical path difference of 2n eff Lo, where n efT is the effective refractive index of the fiber core and L 0 the length of the 
reference cavity. The signals are driven by the tuning of the laser and detected at the photo-detector Dj. They are used 
to trigger the sampling of signal at D 2 , which is the output of another in-fiber interferometer formed with the fiber 
coupler C 2 , a broadband reflector, and a particular fiber Bragg grating at a distance of L,. If there is a series of low 
reflectivity Bragg gratings written at the same wavelength on a single fiber at different locations, the reflected signals 
from each grating are superimposed and detected at D 2 . The detected signals are further processed (fast-Fourier- 
transformed) to obtain the spatial spectrum of all gratings, which displays the physical profiles of the gratings at different 
locations. The spatial spectrum of a particular grating can then be windowed for investigating an individual grating. 
Fig. 2 shows the spatial spectrum of two FBGs in a fiber with multiple gratings. The gratings have a physical length of 
about 5 mm and a separation distance of 10 cm. The signals from the demodulation system are further processed for the 
strain (thermally induced strain, in this case) of each point within the physical length of a grating. Averaging over a 
certain number of points within the grating length is taken for the strain of the particular grating. 



broadband FBGS 
reflector 


Fig. 1. Schematic diagram of a frequency domain demodulation system. Items C 0 , Ci, and C 2 are fiber couplers. “X” 
indicates that the unused port is tenninated. Items D, and D 2 are detectors and FRMs are Faraday rotation mirrors. 
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Fig. 2. The spatial spectrum of two FBGs with a physical length of about 5 mm and a separation distance of 10 cm. 

In this study, an optical fiber with 32 FBGs was bonded to an investigated composite specimen with Kapton® tape, as 
shown in Fig 3 (a). This bonding technique allowed the optical fiber to be bonded and taken off the sample quickly and 




neatly. At the same time it also allowed the fiber to sense the thermal strain for a temperature change up to 50° C 
without imminent slipping. 

The specimen used for testing the viability of this measurement technique was a composite panel with 20 delaminations 
at specified depths. The 10 ply quasi-isotropic composite panel with a lay-up of [0,45, 90, -45, 0/,45, 90, -45/, 0,90/] was 
31.75 x 31.75 centimeters and 0.19 centimeters thick (Fig. 3(a)). The delamination defect areas were squares with sizes 
of 3. 8x3. 8, 2. 5x2. 5, 1.9x1. 9, and 1.3x1. 3 square centimeters. The defects were buried at depths of 10, 20, 30, 40, and 50 
percent of the total thickness. A schematic of the defect layout and FBG sensor numbers is shown in Fig. 3 (b). The 
optical fiber bonded to the specimen had FBGs with a separation (center-to-center) distance of 5.08 centimeters, which 
was the same as the separation distance of the squares in the same row. One FBG was placed at the center of each 
square. Five FBGs numbered 14-19 were outside the squares for comparison. 



Fig. 3. (a) A single optical fiber with multiple Bragg grating sensors bonded onto the surface of a 10-ply-composite, (b) The 
fabricated delaminations of various sizes arranged in columns between plies 1-2 (A), 2-3 (B), 3-4 (C), 4-5 (D), 5-6 (E). 


A quartz lamp was used to heat the front and the back surface of the specimen with a single or multiple heating cycles. 
The cyclic heating had typically a 2 second period with a 50% duty-cycle, and a total of about 20 second duration. The 
grating data acquisition was performed during and following the application of the heat flux to the surface at a rapidity 
rate of 100 Hz. An IR camera was also used to record the temperatures of the surface. The experimental setup is shown 
in Fig. 4. 


4. RESULTS AND DISCUSSION 

Ideally, it is desirable to acquired data from all of the fiber optic sensors at one time. At the higher acquisition rate used 
for these measurements (100 Hz) this was not possible. Instead it was only possible to collect data along 5 adjacent 
gratings. For this reason the subsequent discussion will be limited to data acquired on gratings 8 to 12. These gratings 
are over delaminations at 5 different depths in the composite, with the depth of the delamination increasing 
monotonically from gratings 8 to 12. The time responses during a 1/2.048 Hz heating cycling of the composites as 
measured by gratings 8 to 12 are shown in Fig. 5. As can be seen from this figure, there is no clear trend in the amplitude 
of the thermal responses. This is in part a result of the gratings not being calibrated to enable an absolute temperature 
measurement. 
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Fig. 4. The experimental setup of the fiber optic thermal health monitoring system with a quartz lamp as the heating source 
and a conventional IR camera. 
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Fig. 5. The temperature measurements for 5 different gratings on the composite specimen as result of cycling a quartz lamp 
with a 1/2.048 Flz square wave. Each grating is over a different depth delamination in the composite. 






Fig 6. Linear least squares fit of the thermal response of the composite measured at FBG 8. 


In section 2.2, it was noted that the phase of the response was a function of the thickness and thermal diffusivity of a 
layer and the frequency of the excitation; therefore it is of interest to find the phase of the response for each of the 
gratings. The phase angle has an additional advantage, since it is not dependent on the temperature amplitude, different 
gratings can be compared without calibration 10 . The cycling of power to the lamps does not result in a negative heat flux 
required for a purely cyclic heating of the specimen. Rather, a cyclic heating is superimposed on a static term. The static 
term results from an increase in the temperature of the panel during the test that needs to be accounted for to accurately 
determine the phase. To simultaneously account for the increase in temperature of the panel and calculate the phase, the 
data is fit with 


T(t) = Aj+ A 2 t+ A 3 f 2 + A 4 cos(cot) + A 5 sin(a>t) (8) 

where to is the known angular frequency of the excitation 27i/2.048 Hz. A typical fit of the measured temperature is 
shown in figure 6. The phase angle is calculated from the arctan(A 5 /A 4 ). The measured phase angle for the different 
depth of delaminations is shown in Fig. 7. As can be seen from the figure, there is a large phase angle for delaminations 
closer to the surface. 

For comparison with the measured phase angles, the single layer solution (Eq. (4)), the convolved thermal response (Eq. 
(7)) and a three-dimensional finite element routine were used to calculate the phase angle for different depth 
delamination. The cyclic heating was not synchronized with the data acquisition; therefore the average phase shift for the 
different comparisons was set equal to the measured values. Common literature values of thermal properties of the 
composite were used. The thermal conductivity perpendicular to the fibers was assumed to be 0.97 W irf 1 K/ and 4.85 
W m^K" 1 along the fiber. The specific heat and density were assumed to be 1270 Jkg^K' 1 and 1490 kgm 1 respectively. 
A heat transfer coefficient of 8.063 W m' 2 K'' for convection loss at the front and back surface was used, based on the 


exponential decay in temperature at later times as measured by a fiber optic sensor at the center of the panel. For 
delaminations, 1 .2 cm wide squared delaminations were placed at depths that corresponded to the interface between plies 
at positions corresponding to the squares in the schematic of the composite shown in Fig. 3(b). The smaller 
delaminations have the largest change in phase relative to the one-dimensional model due to in-plane heat flow. 
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Fig 7. Phase measurement for different depths of measurements. The measured phase angles are compared to different 
model results from single-layer soultion, one-dimensional two-layer model, and three-dimensional finite element 
model. 


The comparison of the model results to the measured results is shown in Fig. 7. As can be seen from the figure, the phase 
depends on the contact resistance (gap thickness/thermal conductivity of air) as well as the depth of the delamination. 
The best agreement between models and the measured phase angles is found assuming the air gap in the delamination to 
be 30 pm. Both the one-dimensional model (Eq. (7)) and the three-dimensional finite element model are in reasonably 
good agreement when assuming a 30 pm air gap. This suggests there is relatively little in-plane heat flow in this time 
frame. 

The single layer model shows the same trend as the measured phase angles, however it tends to overestimate the effect 
of depth on the phase angle. This appears to be a result of the finite contact resistance of the delamination significantly 
changing the phase angle. This is in agreement with both the one-dimensional and three-dimensional models for an air 
gap of 0.19 mm, which is a relatively large contact resistance. The estimated phase angles for these models are in good 
agreement with each other and the single layer model. It is clear from these measurements, an estimation of the contract 
resistance is required for accurate depth measurements from the phase measurement. 


5. CONCLUSION 


We have developed a new technique using distibuted FBG sensors for thermographic detection of flaws in materials and 
structures. Individual fibers with multiple FBGs employed as temperature sensors can be bonded to the surfaces of 
structures or embedded in the structures. By applying a cyclic thermal heat flux to the surface of the investigated 
structures, the individual Bragg grating sensors successfully measured the temporal and spatial temperature variations on 
the surface. These thermal responses were consistent with both one-dimensional and two-dimensional models of the 
thermal response for a delaminated composite. The measured phase shifts for delaminations of different depths are in 
good agreement with the model predictions assuming the delamination were 30pm air gaps. Future efforts will focus on 
developing the technique with a faster detection system, detecting more FBGs at one time, and assessing its potential for 
performing thermal health monitoring of aerospace structures and materials. 
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